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We study the geometric phase of a two-level system under the presence of a structured environ¬ 
ment, particularly analysing its correction with the ohmicity parameter s and the onset of non- 
Markovianity. We firstly examine the system coupled to a set of harmonic oscillators and studied 
the decoherence factor as function of the environment’s ohmicity parameter. Secondly, we propose 
the two-level system coupled to a non-equilibrium environment, and show that these environments 
display non-Markovian effects for all values of the ohmicity parameter. The geometric phase of the 
two-level system is therefore computed under the presence of both types of environment. The cor¬ 
rection to the unitary geometric phase is analysed in both, Markovian and non-Markovian regimes. 

Under Markovian environments, the correction induced on the system’s phase is mainly ruled by 
the coupling constant between the system and the environment, while in the non-Markovian regime, 
memory effects seem to trigger a significant correction to the unitary geometric phase. The result is 
significant to the quantum information processing based on the geometric phase in quantum open 
systems. 


I. INTRODUCTION 

A system can retain the information of its motion when 
it undergoes a cyclic evolution, in the form of a geomet¬ 
ric phase, which was hrst put forward by Pancharatnam 
in optics [T] and later studied explicitly by Berry in a 
general quantal system [5]. Since then, great progress 
has been achieved in this field. For example, the applica¬ 
tion of the geometric phase has been proposed in many 
fields, such as the geometric quantum computation. In 
this line of work, many physical systems have been inves¬ 
tigated to realise geometric quantum computation, such 
as NMR (Nuclear Magnetic Resonance) [3], Josephson 
junction [I], Ion trap [5] and semiconductor quantum 
dots [5]. The quantum computation scheme for the geo¬ 
metric phase has been proposed based on the Abelian or 
non-Abelian geometric concepts, and the geometric phase 
has been shown to be robust against faults in the pres¬ 
ence of some kind of external noise due to the geometric 
nature of Berry phase [Mq]. Then, for isolated quan¬ 
tum systems, the geometric phase is theoretically per¬ 
fectly understood and experimentally verihed. However, 
it has been shown that the interactions play an impor¬ 
tant role for the realisation of some specific operations. 
As the gates operate slowly compared to the dynami¬ 
cal time scale, they become vulnerable to open system 
effects and parameters’ fluctuations that may lead to a 
loss of coherence. Consequently, the study of the geomet¬ 
ric phase was soon extended to open quantum systems. 
Following this idea, many authors have analysed the cor¬ 
rection to the geometric phase under the influence of an 
external thermal or non-equilibrium environments, using 
different approaches (see [IIHIS])- In all cases, the purely 
dephasing model considered was a spin-1/2 particle cou¬ 
pled to the environment’s degrees of freedom through a 
CTjj-coupling. The interest on the geometric phase in open 
systems has also been extended to some experimental se¬ 


tups m- Lately, it has also been observed in a vari¬ 
ety of superconducting systems [HUH!, and it has been 
shown the importance of quantifying decoherence when 
geometric operations are carried out in the presence of 
low-frequency noise. All real experiments generally im¬ 
ply the presence of an external environment which in¬ 
duces noise and dissipation on the subsystem depending 
on the strength of the coupling among them. Memory 
effects are also considered a relevant source of noise that 
can affect the dynamics of the system of interest. Thus, a 
detailed mechanism of decoherence due to external noise 
sources is still required in order to overcome the effect of 
a destructive decoherence in measurements of geometric 
phase. 

Within a microscopic approach, quantum Markovian 
master equations are usually obtained by means of 
the Born-Markov approximation, which assumes a weak 
system-environment coupling and several further, mostly 
rather drastic approximations m- However, these ap¬ 
proximations are not applicable in many processes oc¬ 
curring in nature, such as strong environment couplings, 
structured and finite reservoirs, low temperatures and 
large initial system-environment correlations [2Ij . In 
the case of any substantial deviation from the dynam¬ 
ics of a quantum Markov process, one often speaks of 
a non-Markovian process, implying that the dynam¬ 
ics is governed by significant memory effects. Unfor¬ 
tunately, a consistent general theory of quantum non- 
Markovianity does not exist and even the very dehni- 
tion of non-Markovianity is a highly up-to-date issue. 
Very recently important steps towards a general theory 
of non-Markovianity have been made. There has been 
a great effort to rigorously define the border between 
Markovian and non-Markovian quantum dynamics and 
to develop quantitative measures for the degree of mem¬ 
ory effects in open systems [221 - 124] . The quantification 
of non-Markovianity is justified by the fact that there is 
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increasing evidence of its crucial role as a resource for 
quantum technologies [24l [25]. Non-Markovianity evo¬ 
lution is often characterised by decoherence phenomena 
and information trapping, hence leading to longer co¬ 
herence times in comparison to the Markovian case. As 
reservoir engineering techniques become experimentally 
feasible, it is crucial to establish relations between the 
occurrence of non-Markovianity and the form of the en¬ 
vironmental spectrum. Recently, there has been many 
studies on systems whose reduced dynamics are charac¬ 
terised by memory effects. Non-Markovian features play 
an important role in systems where the frequency spec¬ 
trum of the environment is structured. In Ref. |2hj . au¬ 
thors established a connection between the general form 
of the spectrum and the memory effects in the reduced 
system dynamics. 

In this framework, this study has been two-folded mo¬ 
tivated: on the one side, the need to have a better under¬ 
standing of the geometric phase in open quantum systems 
in order to achieve fault tolerance quantum computation. 
On the other side, the need to understand memory ef¬ 
fects as a source of noise for the system of interest. In 
this manuscript, we shall study the correction to the geo¬ 
metric phase of a two-level system coupled to an external 
environment. We shall use the spin-boson model since it 
has the advantage of having an analytical solution which 
can be useful to have a better insight into the onset of 
non-Markovian effects. By the use of the established re¬ 
lation between the form of the reservoir spectrum and the 
onset of non-Markovianity, we shall study the geometric 
phase of the open system as a function of the “ohmic- 
ity” of the reservoir, which allows the description of sub- 
Ohmic, Ohmic and super-Ohmic spectrum. The paper is 
organised as follows: in Section II, we briefly described 
the model and present both types of environments to be 
analysed: thermal equilibrium and non-equilibrium ones. 
We study the diffusion coefficients and the corresponding 
decoherence factors. Following Ref.|5S| we study the on¬ 
set of non-Markovianity in the environments considered. 
In Section III, we compute the geometric phase acquired 
by the system for different forms of the reservoir spec¬ 
trum, as function of the ohmicity parameter and study 
how the memory effects of the environment affect the ge¬ 
ometric phase. Finally, in Section IV, we shall make our 
final remarks. 


II. THE MODEL 

The spin-boson model is studied in a variety of fields, 
such as condensed matter physics, quantum optics, phys¬ 
ical chemistry and quantum information science [27] in 
order to describe non-unitary effects induced in quantum 
systems due to a coupling with an external environment. 
For a quantum system, the influence of the surroundings 
plays a role at a fundamental level. When the environ¬ 
ment is taken into consideration, the system dynamics 
can no longer be described in terms of pure quantum 


states and unitary evolution. From a practical point 
of view, all real systems interact with an environment, 
which means that we expect their quantum evolution to 
be plagued by non-unitary effects, namely dissipation and 
decoherence. Most theoretical investigations of how the 
system is affected by the presence of an environment have 
been done using a thermal reservoir, usually assuming 
Markovian statistical properties and defined bath corre¬ 
lations [23 123 (there are also works on non-Markovian 
models as, just for example, [30]). 

In the following, we shall consider a paradigmatic spin- 
boson model consisting of a two-level system coupled to 
an external environment. We shall consider two differ¬ 
ent types of environments and see the non-Markovian 
effects: thermal equilibrium and non-equilibrium envi¬ 
ronment. In both cases, we shall compute the diffusion 
and decoherence factor derived by the definition of the 
corresponding bath correlations and the spectral density 
/(w), which can be defined for a general environment as: 

H^) = 70 ^^ exp(-w/A). (1) 

In Eq.Q, A is the reservoir frequency cutoff and 70 is 
the coupling constant (has different units for the different 
environment considered). By changing the value of the 
s-parameter one goes from sub-Ohmic reservoirs (s < 1) 
to Ohmic (s = 1) and super-Ohmic (s > 1) reservoirs, re¬ 
spectively. The Ohmic spectrum is generally used to de¬ 
scribe charged interstitials in metals. The supra-Ohmic 
environment commonly describes the effect of the inter¬ 
action between a charged particle and its electromagnetic 
field. The sub-ohmic environment is used to model the 
type of noise occurring in solid state devices due to low 
frequency modes, similar to the “1/f” noise in Josephson 
junctions. Thus, the description of the spectral density 
function in terms of a continuous parameter s allows to 
simulate paradigmatic models of open quantum systems 
by the variation of s. We stress that such engineering of 
the ohmicity of the spectrum is possible when simulating 
the dephasing model in trapped ultracold atoms |31j . 


A. Thermal equilibrium environments and the 
onset of non-Markovianity 

We shall start by studying the decoherence effects in¬ 
duced by a thermal equilibrium environment, generally 
modelled by a set of harmonic oscillators. The interaction 
between the two-state system and the environment is en¬ 
tirely represented by a Hamiltonian in which the coupling 
is only through az- In this particular case, [cr^, Hint] = 0 
and the corresponding master equation for the reduced 
density matrix, is much simplified, with neither frequency 
renormalisation nor dissipation effects. The Hamiltonian 
of the complete system is written as: 
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Hsb = Afc(gfcafc + fffcgfc) + hujkalak- 

k k 

( 2 ) 

The interaction hamiltonian is only proportional to Uz 
which means that the populations of the eigenstates re¬ 
main constant while the off-diagonal terms of the reduced 
density matrix decay due to the existence of the environ¬ 
ment, as 


where J- is called the decoherence factor defined as (see 
Ref. [13] for details) 


poo 

T{t) = 2 / da;/(w)coth 

Jo 

= 2 / dsD{s), 

Jo 


w \ (1 — cos(a;t)) 


2kBT 


( 4 ) 


where D{s) is the diffusion coefficient, present in the 
master equation. It is clear that not only depends 
on the temperature of the bath but also on the spectral 
density, particularly on the ohmicity parameter s. 


In Ref. [331 , authors have shown that there exists a link 
between the onset of non-Markovianity and the form of 
the reservoir spectrum for thermal equilibrium environ¬ 
ments as the ones considered here. They have stated 
that the crossover is signalled by the onset of periods 
during which the diffusion rate is negative. A common 
feature of all non-Markovianity measures is that they are 
based on the non-monotonic evolution in time of cer¬ 
tain coefficients which signal the information back-flow 
from the environment back to the main system. Without 
memory effects, decoherence rates (as a quantum chan¬ 
nel) decreases in time monotonically. However, environ¬ 
ment memory effects may produce a non-monotonic be¬ 
haviour of the quantum channel capacity. This is related 
to the convex to non-convex changes of the spectrum. 
The condition on the non-convexity of the environmen¬ 
tal spectrum is a necessary and sufficient condition for 
non-Markovianity at all temperatures. This has been set 
as s > 2 for zero-T environments. Hence, memory ef¬ 
fects leading to information back flow and recoherence 
occur only if the reservoir spectrum is super-Ohmic with 
s > 2. This means that even if the reduced dynam¬ 
ics is exact, and hence no Markovian approximation has 
been performed, the time evolution of the qubit does not 
present any memory effects for Ohmic and sub-Ohmic 
spectra. This argument has been derived based on the 
non-monotonic behaviour of the decoherence factor J^(t), 
and basically it is equivalent to the fact that the diffusion 
coefficient D{s) becomes negative. 

In Fig{TJ we present a contour plot for the diffusion 
coefficient for different equilibrium environments at zero 
temperature, to show the onset of non-Markovianity as 
derived in [33] • We can see the regions where the diffusion 


coefficient has a negative value, which are related to the 
onset of non-Markovian effects. For example, in Figj^ 
we see that for Ohmic thermal environments, i.e. s = 1 
in the vertical axis, we always have positive values for 
the diffusion coefficient. This means that the Ohmic en¬ 
vironment is always Markovian. For other environment, 
such as s = 4 in the vertical axis, we can find regions 
where the diffusion coefficient gets positive values and 
areas where the value is negative. This means that the 
memory effects become important after some time evo¬ 
lution. 



FIG. 1. (Color online) Density plot of the diffusion coefficient 
of thermal environments. The horizontal axis represents the 
time evolution as Q.t and the vertical axis is the ohmicity s. 
The darkest areas are negative values of the diffusion coef¬ 
ficient which has been shown to be related to the onset of 
non-Markovian effects as stated in [^. Parameters used: 
7 o =0.1 and A = lOH. 


Using Eqs. Q and Q, the exact decoherence factor 
as function of time and the ohmicity parameter s, can be 
written as: 


■^(t) = 47o^ fl-(l+A2t2) 4 
s — 1 L 

X (cos[s arctan(At)] -h At sin[s arctan(At)])], (5) 


where r[a:] is the Gamma function. In Figj^ we present 
the evolution of the decoherence factor as time evolves. 
In Ohmic (s = 1) thermal environments, the decoherence 
factor is known to be a monotonic decreasing function in 
time (blue line in Figj^. However, when the memory 
effects of the environment become relevant (for example 
s = 4), the decoherence function has not a monotonic be¬ 
haviour, leading to recoherence phenomena (red dashed 
line in Figj^. This is related to the negative regions of 
the diffusion coefficient as can be seen in Fig0 
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FIG. 2. (Color online) Decoherence factor as function of the 
time for different values of the ohmicity parameter s. It can 
be seen that areas with negative diffusion coefficient in Fig[^ 
present a non-monotonic behaviour in the decoherence factor 
leading to recoherence effects. The blue line corresponds to 
s = 1 and the red dashed one corresponds to s = 4. Parame¬ 
ters used: 70 = 0.1 and A = 1011. 


B. Non-equilibrium environments and the onset of 
non-Markovianity 

In this subsection, we shall deal with non-equilibrium 
environments, represented by random perturbations with 
non-stationary statistics. The motivation to study these 
type of environments is twofold: on one side, the diffu¬ 
sion coefficients computed for different environment spec¬ 
tra have negative parts for all values of the ohmicity 
parameter. This means they can be considered as non- 
Markovian environments for all values of the parameter s 
and, contrary to the environments studied above, there is 
not a critical value of s which determines the onset of non- 
Markovianity. On the other side, these non-equilibrium 
baths may represent a proposal for engineering reservoirs 
in a manner reminiscent of a coherent control experiment 
using shaped pulses [32]. In this model, the control pa¬ 
rameter A, is derived not from a laser pulse, but rather 
from well-defined phase relations between the modes of 
the bath. 

The modelling of these non-equilibrium environments 
implies that the two-level quantum system presents an 
energy gap E 2 {t) — Ei{t) = huj{t) which fluctuates due 
to the influence of the environment, where Ej{t), with 
j = 1,2 is the instantaneous energy of state j as per¬ 
turbed by the surroundings. Following the idea pro¬ 
posed in |33|, the bath is represented by a random func¬ 
tion of time corresponding to the transition frequency 
of the two-state quantum system. In contrast to the 
usual treatment, the statistical properties of this random 
function are non-stationary, corresponding physically to, 
for example, impulsively excited phonons of the envi¬ 
ronment with initial phases that are not random, but 
which have defined values at t = 0. Due to these as¬ 
sumptions, this environment is not at thermal equilib¬ 
rium. The time-dependent frequency is written in the 


form uj{t) = il + 6oj{t). 



FIG. 3. (Golor online) Density plot of the diffusion coefficient 
for the non-equilibrium environments. The horizontal axis 
represents the time evolution as fit and the vertical axis is 
the ohmicity s. The darkest areas are negative values of the 
diffusion coefficient. Parameters used: A = lOD , 70 = 0.1, 
DA = 0.3 and d = 2D. 

Following this approach, it is easy to check that the 
off-diagonal terms of the reduced density matrix can be 
written as 

= Prc,(0). (6) 

Here, we denote with < ... > the non-equilibrium av¬ 
erage over the non-stationary random bath and E(t) is 
defined as the decoherence factor for the non-equilibrium 
environments, which reads, 

T{t) = exp{— (70 exp(—4c?f)(—1 -I-exp(2(it) 

X r[l + s]((l + 4(t-A)2A2)^^ 

X cos((l + s) arctan(2A(t — A))) + cosh{2dt) 

+ sinh(2d<)))}, (7) 

where d and A are parameters of the environment model. 
By knowing the decoherence factor we can have a better 
insight into the decoherence process induced in the sys¬ 
tem by the presence of a non-equilibrium environment. 

In Figj^ we show a contour plot for the diffusion co¬ 
efficients for different forms of the ohmicity parameter s 
in the vertical axis and different values of fit in the hor¬ 
izontal one. We can note that there is always a period 
when the coefficient becomes negative, even for s = 1. 
This means that memory effects are always present for 
this type of environments and the decoherence factors 
are always non-monotonic decreasing functions, allowing 
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FIG. 4. (Color online) Decoherence factor for different values 
of the ohmicity parameter s as function of time. We can see 
that the dip becomes more relevant and drastic for bigger 
values of s. The red dashed line is for s = 1, the black solid 
line is for s = 2 and the blue dotted line for s = 3. Parameters 
used: A = 1012 ,70 = 0.1, DA = 0.3 and d = 2Q,. 

recoherence effects on the system of interest. We can 
consider these environments to be non-Markovian for all 
values of s. We can also note a clear hierarchy in the 
coefficients: the bigger the value of s, the more nega¬ 
tive is the diffusion coefficient and the more decoherence 
(and recoherence) induced by the environment. This fact 
can be observed in Fig|^ where we show the decoher¬ 
ence factor for different non-equilibrium environments. 
It can be seen that even though we have s = 1, the 
behaviour is non-monotonic, quite different to the deco¬ 
herence factor of a thermal equilibrium environment. In 
Figj^we present the decoherence factor for different val¬ 
ues of the ohmicity parameter s as function of time. The 
red dashed line corresponds to s = 1, the black solid line 
is for s = 2 and the blue dotted line shows the case s = 3. 
In all cases, these non-equilibrium environments present 
a “dip”, which makes the behaviour of the decoherence 
factor non-monotonic. The strength and location of the 
“dip” is determined by the other parameters of the model 
A and d. 


where efc(t) are the eigenvalues and Idffe) the eigenstates 
of the reduced density matrix (obtained after tracing 
over the reservoir degrees of freedom). In the last defini¬ 
tion, T denotes a time after the total system completes a 
cyclic evolution when it is isolated from the environment. 
Taking into account the effect of the environment, the 
system no longer undergoes a cyclic evolution. However, 
we shall consider a quasi cyclic path V t e [0, r], with 
T = 27r/r2 (H is the system’s characteristic frequency). 
When the system is open, the original GP, i.e. the one 
that would have been obtained if the system had been 
closed (j)u, is modified. This means, in a general case, that 
the phase can be interpreted as (I>g = (l>u + d4>, where Scj) 
depends on the kind of environment coupled to the main 
system [T3HIS1II71IM1I35]. 

We want to compute the GP acquired by the system 
for different forms of the reservoir spectrum, as function 
of the ohmicity parameter, which allows the description 
of sub-Ohmic, Ohmic and super-Ohmic spectra. Partic¬ 
ularly, we want to see how the unitary geometric phase 
for a two-level system (fiu = 7J'(1 — cos 9) is corrected as a 
function of the ohmicity, i.e. the parameter s of the spec¬ 
tral density /(w). Assuming an initial quantum state of 
the system as 

IV’(O) >= cos(6</2)|0 >-I- sin(0/2)|l >, (9) 

its evolution at a later time t, is 

IV'(I) >= cos(0+(t))|O > -I- sin(0+(t))|l >, (10) 

where cos(0+(t)) (and sin(0-|_(t))) encodes diffusion in¬ 
duced on the subsystem due to the presence of the en¬ 
vironment. As explained in Ref. m, the GP is obtained 
by computing eigenvectors and eigenvalues of the reduced 
density matrix derived by using the state vector \tp{t) > 
and using Eq.([^. 

We shall consider the thermal equilibrium and non¬ 
equilibrium environments considered above. With the 
decoherence factors and the reduced density matrix com¬ 
puted analytically for each case, we can compute the GP 
and see how it is affected by the memory effects present 
in the different environments considered. 


III. CORRECTION TO THE GEOMETRIC 
PHASE 


A. Thermal equilibrium environments 


In order to compute the geometric phase (GP) and 
note how it is corrected by the presence of the environ¬ 
ment, we shall briefly review the way the GP can be 
computed for a system under the influence of external 
conditions such as an external bath. In Ref. m , a quan¬ 
tum kinematic approach was proposed and the GP for a 
mixed state under non-unitary evolution has been defined 
as 

4 >g = arglXl V£^fe(0)£fc(T)(4'fc(0)|T'fe(T)) 

k 

( 8 ) 


By considering the thermal equilibrium environments 
in the zero-T limit and using the corresponding decoher¬ 
ence factor (Eq. (§), we can find an expression for the 
correction to the GP in terms of a series expression in 
the coupling constant 


= 4’u+ 7o sin^ 0COS0 r(s — 2)|27r(s — 2) 


dTT^A^N - ! 


27rAN 


47r"-A^\-2r, ( (I-kLVW 

+ ( 1 +-fi^) (»“«“(-^)) 


D2 

att^a n ( 

+ (—-X>“"( 




s arctan 


(“))]}■ 


( 11 ) 
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As it is expected, in the limit s —> 1 the 
correction to the phase approaches 5(j) « 

47 r 7 o sin^ 0 cos 6 * (—1 + log[27rA/f2]). In the limit s —>■ 3, 
the correction is given by 5(j) « dvryo sin^ 0 cos 0. Both 
expressions agree with the corrections to the Ohmic 
(s = 1) and supra-Ohmic (s = 3) geometric phases in the 
zero temperature limit found in Ref. m- As this result 
is derived for small values of the coupling constant 70 , in 
the following we shall compute the exact GP numerically. 

In Fig[^ we show the correction induced by different 
environmental types on the geometric phase (normalised 
by the value of the unitary geometric phase (pu) as a 
function of the ohmicity for different values of 70 . In this 
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FIG. 5. (Color online) Correction to the geometric phase (in 
units of the unitary phase) as a function of the ohmicity pa¬ 
rameter s. Red cross-solid line is for 70 = 0.001, blue asterisk¬ 
line is for 7 o = 0.005, magenta square-line is for 70 = 0.01 and 
black circle-line for 70 = 0.03 for zero T environments. Pa¬ 
rameters used: 0 — tt/S, A = lOfl. 

figure we can note that the geometric phase is very much 
destroyed when s —>■ — 1 . In that case, we are considering 
the effect of a noise, similar to 1//, which is very harm¬ 
ful. This type of noise can be considered a sub-Ohmic 
environment due to the fact that low frequencies are pre¬ 
dominant. The correction to the phase grows as the effect 
of low frequency modes of the environment becomes more 
relevant. On the contrary, as s starts increasing, we ob¬ 
tain the correction to the phase similar to the one of an 
Ohmic environment (s = 1 ) at zero-T [T^. As can be 
expected, decoherence induced by this type of environ¬ 
ment is low for a very weak coupling (less than 10% for 
smaller values of 70). However, it becomes significant for 
bigger values of 70, for example for 70 = 0.03 in Figj^ 
the correction is bigger than 20%. This agrees with the 
results in Ref. [ 13 ] : more decoherence induced on the sys¬ 
tem, implies a bigger correction to the unitary geometric 
phase. Near the Ohmic region of the parameters set, the 
correction depends mainly on 70, as expected, being neg¬ 
ligible when the small value of 70 is such that is not able 
to destroy coherences in the system. For 1 < s < 2 the 


correction to the phase is more insensitive to the ohmicity 
value, and it is evident that there is a change in behaviour 
for s > 2 (it starts increasing), in agreement with the on¬ 
set of non-Markovianity. This shows that non-Markovian 
environments induce a bigger correction to the unitary 
geometric phase. As expected, this fact is enlarged for 
bigger couplings between system and environment, i.e. 
values of 70 , since these imply a bigger decoherence ef¬ 
fect. These corrections have not been studied previously, 
for example in Ref. [13]. In conclusion, corrections to the 
phase are double. Firstly, we have the common hierar¬ 
chy in the induced correction ruled by the coupling to 
the environment: the bigger correction occurs with the 
bigger value of 70 . Secondly, it is possible to see that 
the correction to the phase becomes more relevant when 
the ohmicity parameter surpasses its critical value. It 
becomes evident that the more non-Markovian the envi¬ 
ronment, the biggest correction to the geometric phase 
for the same coupling to the environment. This explains 
the behaviour observed in Figj^ 

In Figj^the normalised correction is plotted for differ¬ 
ent environments at zero temperature. Red asterisks rep¬ 
resent an Ohmic environment (s = 1), while blue squares 
are for a supra-Ohmic environment of s = 2 and black 
circles one of s = 2.5. The triangles represent s = 3. 
We can see that for s > 2, the correction to the geomet¬ 
ric phase is bigger and its behaviour is more drastic, in 
agreement with the onset of non-Markovianity expected 
to be for environments at zero temperature. 
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FIG. 6 . (Color online) The correction to the GP as function of 
the angle that defines the state vector in the Bloch sphere (6) 
for Ohmic and supra-Ohmic environments in the zero T limit. 
Red asterisks represent the Ohmic environment (s = 1), while 
blue squares are for a supra-Ohmic environment of s = 2 and 
black circles one of s = 2.5. The triangles represent s = 3. 
Parameters used: 70 = 0.01 and A = 1017. 
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B. Non-equilibrium environments 


Finally, we compute the correction to the GP for the 
non-equilibrium environments which have memory effects 
for all values of s. Therefore, we shall consider them as 
non-Markovian environments. 

In Figj^ we show the behaviour of the normalised cor¬ 
rection to the unitary geometric phase as a function of 
the ohmicity s. In this figure we can see that for smaller 
values of s and particular values of A and d parameters 
(which determine a small dip in the decoherence factor 
coefficient as shown in Figj^, the non-equilibrium envi¬ 
ronment does not have a big influence on the geometric 
phase, and the open geometric phase coincides with the 
unitary one. In Fig. we present different sets of values 
for the parameters of the environment’s model. It is easy 
to note that the correction is mainly ruled by the value of 
the coupling constant: see, for example, the blue dotted 
line and the blue dots. Both lines correspond to the same 
value of 7 o = 0.5 but different frequency cutoff A. The 
same occurs with the red line and the red asterisks with 
7 o = 0.1. The magenta triangles share the same value of 
7 o but the values of A and d are interchanged with respect 
to the red line and asterisks. Finally, the black lines in¬ 
dicates smaller values of 70 . When s increases (even for 
the set of parameters 70 , A and d, for for which deco¬ 
herence is negligible) the effect of the memory effects on 
the geometric phase of the non-equilibrium environment 
are stronger and the correction to the geometric phase 
results bigger. After this, the non-Markovian correction 
to the phase presents an abrupt slope, leading to bigger 
corrections compared to the equilibrium baths considered 
above. By an analytical comparison of the perturbative 
Eqs. (11) and (12), it can be seen that the correction in¬ 


duced by non-equilibrium environments is similar to that 
of the thermal ones for ohmicity parameters — 1 < s < 1 
(for equal coupling constant 70 and environment cutoff 
A). However, as s increases, the correction induced by 
non-equilibrium environments gets bigger in comparison, 
becoming particularly important for s > 3. This confirms 
the importance of the memory effects of the environment 
in the correction of the geometric phase, as already stated 
in Refs. [55] and |35| , where authors studied the memory 
effects of the environment for different two-level systems. 

As in previous Section, we can estimate the correction 
to the unitary geometric phase, in an expansion in pow¬ 
ers of the coupling between system and environment 70 . 
In this case, the geometric phase can be approximated, 
analytically as 


4>G ^ 4>u+7oT^[s + 1] sin^ OcosO , 


( 12 ) 


where we have neglected 0 ( 7 oH/A) terms. This ap¬ 
proximate expression matches very well the low -70 
corrections given in Fig. This expression also gives 
the leading correction to the phase found in Ref. |T 6 | for 
s = 1 and s = 3. On general grounds, we can see that 
the non-equilibrium bath is more harmful and induces a 
bigger correction on the geometric phase for all values of 



FIG. 7. (Color online) Correction to the unitary geometric 
phase for different values of the model’s parameters.The blue 
dotted line is for 70 = 0.5 and A = lOH while the blue square 
dots is 70 = 0.5 and A = 5H. The red asterisks are for 
70 = 0.1 and A = 5H, while the solid red line is 70 = 0.1 and 
A = in and the triangle dotted line 70 = 0.1 and A = IH 
but changing the parameters of A and d. Finally, the dotted 
black line is for 70 = 0.01 and A = lOH. Parameters used: 
HA = 0.5, d = IH. 


0 than the equilibrium environment. 


IV. CONCLUSIONS 

The geometric phase of quantum states could have a 
potential application in holonomic quantum computation 
since the study of spin systems effectively allows us to 
contemplate the design of a solid-state quantum com¬ 
puter. However, decoherence is the main obstacle to over¬ 
come. Furthermore, in most cases of practical interest, 
quantum systems are subjected to many noise sources 
with different amplitudes and correlation times, corre¬ 
sponding de facto to a non-equilibrium environment. 

We have computed the correction of the geometric 
phase under the presence of structured reservoirs. We 
have defined the spectral density as function of the 
ohmicity parameter s. This is advantageous because it 
allows to study sub-Ohmic, Ohmic and supra-Ohmic en¬ 
vironments in the same approach. One could wonder if 
the correction to the geometric phase is due to the non- 
Markovianity of the environment or simply to a stronger 
effective interaction among the system and the environ¬ 
ment. 

Firstly, we have considered the structured reservoirs 
to be composed of a set of harmonic oscillators at zero 
temperature. In the case of thermal equilibrium environ¬ 
ments, we have shown that for small values of the ohmic¬ 
ity s < 2 , the hierarchy on the correction of the unitary 
geometric phase is mainly due to the value of 79 . This 
means: the stronger the coupling to the environment, the 







stronger the correction induced on the geometric phase 
as expected. However, for s > 2, even though the hier¬ 
archy is repeated, we can see that the non-Markovianity 
has also an important role in the correction. We have 
checked this argument by considering the decoherence 
factor for different values of s and observing that the 
correction to the geometric phase is strongly influenced 
by this coefficient. 

Secondly, we have considered the structured environ¬ 
ments to be modelled by a type of non-equilibrium en¬ 
vironments. The non-equilibrium feature is represented 
by a non-stationary random function corresponding to 
the fluctuating transition frequency between two quan¬ 
tum states coupled to the surroundings. We have shown 
that the diffusion coefficients have always negative values 
for same period of time for s > 1 which means a non¬ 
monotonic behaviour of the decoherence factors. This 
allows to consider these non-equilibrium environments as 
non-Markovian ones for all values of s. In this framework. 


we have computed the non-unitary geometric phase for 
the qubit in a quasi-cyclic evolution under the presence 
of these particular non-equilibrium environments, both 
numerically and analytically. When compared the cor¬ 
rection of the geometric phase induced by these envi¬ 
ronments and the thermal ones, memory effects of the 
environments with random noise are more harmful and 
the correction induced is bigger particularly for values of 
s > 3. 

Finally, these kind of environments could become a 
proper experimental setup for the observation of the ge¬ 
ometric phase. Our work on the assessment of the en¬ 
vironmental effect on the geometric phase, especially in 
the non-Markovian regime (either for equilibrium or non¬ 
equilibrium environments) could be of great importance 
in using the geometric phases in two-level systems to im¬ 
plement the quantum gates. The results presented in 
this paper can also provide a clue to observe the GP in 
a two-level system. 
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